Steady state entanglement in open and noisy quantum systems at high temperature 
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We show that quantum mechanical entanglement can prevail in noisy open quantum systems at high temper- 
ature and far from thermodynamical equilibrium, despite the deteriorating effect of decoherence. The system 
consists of a number N of interacting quantum partides, and it can interact and exchange partides with some 
environment. The effect of decoherence is counteracted by a simple mechanism, where system partides are 
randomly reset to some Standard initial state, e.g. by replacing them with partides from the environment. We 
present a master equation that describes this process, which we can solve analytically for small TV. If we vary 
the interaction strength and the reset against decoherence rate, we find a threshold below which the equilibrium 
state is classically correlated, and above which there is a parameter region with genuine entanglement. 

PACS numbers: 03.65.Yz, 03.65.Mn,03.67.-a 



While in quantum computation and communication the rel- 
evance of entanglement is undisputable and a number of ex- 
isting protocols use specific entangled states as resources, the 
situation is not so clear in less controlled, macroscopic sys- 
tems such as sòlids or fluids, consisting of a large number of 
partides, which interact with each other and with the envi- 
ronment. Theoretical studies of so-called spin chains and lat- 
tices show that (bipartite and multipartite) entanglement 
is present in these systems, at least at very low temperatures 
when they are close to their energètic ground state yj. This 
can be qualitatively understood from the fact that the ground 
state of most interacting systems is entangled with respect to 
their interacting constituents. By a continuity argument, then, 
one expects this to be true also for small temperatures above 
zero. The situation is different, however, for gas-type sys- 
tems and systems at higher temperatures, or, more generally, 
systems that are also allowed to exchange partides with the 
environment and that may be far from thermodynamical equi- 
librium. Examples of such systems include man-made devices 
such as the làser but also e.g. bio-molecular processes 

in celis. 

We study open driven quantum systems from the perspec- 
tive of quantum information. The main question is whether 
in interacting systems, where the individual constituents are 
subjected to decoherence, and where the particle numbers 
may fluctuate, entanglement can exist in steady-state. Except 
in situations with very special decoherence mechanisms, or 
in strongly interacting systems at very low temperature, one 
expects that decoherence diminishes and inevitably destroys 
entanglement. Here, we consider a system of N interacting 
quantum partides (qubits, for simplicity) and ask whether en- 
tanglement can exist under the following conditions. (i) The 
interaction between the partides is capable of creating entan- 
glement within the system. (ii) There is an interaction be- 
tween the partides and an environment, which will lead to de- 
coherence on the state space of the system. (iii) The environ- 
ment is not allowed to introduce entanglement, i.e. it is sup- 
posed to interact only locally with the system partides, lead- 
ing to individual decoherence. Furthermore, it cannot "im- 
port" any entanglement into the system (e.g. by exchanging 
two system partides with a pair of entangled partides from 



the environment, nor by some other external control as it is 
e.g. assumed in quantum error correction). (iv) The system 
may exchange partides with the environment, in which case 
the number N may fluctuate around some average value N. 
An example for a system fulfilling (i-iv) is a spin gas 01, 
where partides carrying a spin degree of freedom move freely 
and interact for a short time upon collision. The partides are 
subjected to individual decoherence between two collisional 
events, and also the number of system partides may fluctu- 
ate. We remark that condition (iii) may not be fulfilled in 
certain situations, e.g. for strongly interacting systems cou- 
pled to a thermal bath. In this situation, interactions between 
the compound system and the environment may drive tran- 
sitions to entangled states, allowing for entanglement of the 
thermal state at low temperatures. However, a local decoher- 
ence mechanism as demanded in (iii), which especially takes 
account of gas-type scenarios where system partides move in 
a background gas of other molècules, with which they collide, 
makes it much more difficult to obtain steady-state entangle- 
ment. 

In this paper, we identify a simple mechanism that satisfies 
the criteria given above, and which can still lead to entangle- 
ment in a steady state not necessarily close to the ground state 
or thermal state of the interaction Hamiltonian. The mecha- 
nism corresponds to replacing randomly (and at a given rate r) 
system partides with partides from the environment in some 
Standard, sufficiently pure, single-particle state. The main pur- 
pose of this mechanism is that it locally extracts entropy from 
the system. Thereby, existing entanglement between the to- 
be-replaced particle and the rest of the system is destroyed, but 
the "fresh" particle is ready to become entangled with the sys- 
tem via the interaction Hamiltonian. This simple mechanism 
is equivalent to a situation where each particle is measured 
with the rate r and then reset to a Standard state (ignoring the 
correlations with the other partides). This mechanism, taken 
alone, certainly can not introduce entanglement into the sys- 
tem, as it acts locally on the individual partides. However, in 
combination with the interaction between system partides, it 
is able to create fresh entanglement on an appropriate coarse 
grained time scale, and thus to counteract the effect of de- 
coherence. The mechanism bears some analogy with driven 
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systems such as an atòmic -beam làser, where excited àtoms 
are injected into a làser cavity, providing a resource of energy 
in form of via stimulated emission of photons, thereby main- 
taining the state of the làser field. Similar as in làser theory, 
we find a threshold: If, for a given decoherence and interaction 
strength, we increase the reset rate r, we find a threshold-value 
below which the steady state of the spins is always classically 
correlated, while for larger vàlues of r the system shows quan- 
tum correlations (i.e. entanglement). 

In the language of reservoir theory 0|, such a system is 
coupled to two different reservoirs: one at high temperature 
being responsible for decoherence, and a second at low tem- 
perature providing fresh partides with low entropy. The sys- 
tem is thus far away from thermodynamical equilibrium. Re- 
markably, the proposed local(!) reset mechanism allows one 
to obtain steady state entanglement quite generically, indepen- 
dent of the specific interaction Hamiltonian and of the details 
of system-environment coupling (e.g. temperature of bath), 
and also in case of fluctuating particle numbers. This is in 
contrast to entanglement in thermal states, which is restricted 
to low temperatures (see e.g. 

In the following, we capture the essential features of such 
open and noisy systems in a toy model consisting of only two 
qubits. We give the master equation describing the quantum 
mechanical interaction, individual noise channels, and the re- 
set mechanism, satisfying conditions (i-iii). We solve the mas- 
ter equation and obtain fully analytic expressions for the en- 
tanglement in the system. In a second part, we show that the 
properties of the toy model carry over to larger systems, even 
with fluctuating particle numbers. 

Model. The system consists of N qubits that interact 
with each other according. We describe this interaction on 
a coarse-grained timescale by an effective Hamiltonian H. 
This includes gas-type scenarios with short time interac- 
tions, as well as continuous coupling between system com- 
ponents. Additionally, the partides are subject to decoher- 
ence described by a noise channel and to a reset channel, 
modelling the system's interaction with the environment. We 
describe the noise channel £ no ise by the Lindblad opera- 
tor @: C noise p = Ett -1(1 - s)[afa®p + pafa® - 
2a® pa®] - §s[a®afp + pa®af - 2af pa®] - 

2C ^ B [p — a®pa®], where a± = (a x ± ia y )/2 and the 
crs are Pauli operators. Parameters B and C give the de- 
cay rate of inversion ( ) an d polarization (a±), and s = 
( 2* )<x> e 1] depends on temperature T, where s = 1/2 
corresponds to T = oo. This form of the noise channel is 
based on certain approximations, e.g. Markov approxima- 
tion. Note, however, that this is not an essential assumption 
in the present context, but it simplifies the analytic treatment. 
The reset channel corresponds to an additional term speci- 
fied by Lioville operator C rese t, which we model as follows: 
CresetP = r Ç i=1 2 ( \Xi\i (Xi I tr iP ~ p)- lt means that with 
some probability rot particle i, i = 1...N, is reset during the 
time interval St to some specific state \xí)í, while the other 
qubits are left in the state tiip. The total master equation is 



then given by 

N 

P = -i[H,p] + CnoiseP + ^2r(\xi)i(xi\trip - p). (1) 

i=l 

We have checked that this equation is of Lindblad form 01 
and hence it generates a time evolution of p described by 
a completely positive map. For N = 2, we have analyti- 
cally solved the equation for Heisenberg and Ising Hamilto- 

nians, Hh = ga^ ■ a^ 2 \ Hj = gai az, where g > 
is the coupling strength. In the following we will illustrate 
the essential features for the Ising Hamiltonian Hj, and de- 
coherence by a purely dephasing channel. This amounts to 
setting B = (then s is arbitrary) and C = 7, leading to 

Cnoise = l/2J2i=i,2Í a ^ P a ^ _ p)- Fimilly, we consider 
the case \xí)í = |+)i> where (7 X |+) = |+) is an eigenstate of 
a x . The choice of the reset state \xí)í depends on the Hamil- 
tonian H, which should be able to create entanglement from 
the resulting product state. The Ising Hamiltonian could e.g. 
not create any entanglement if \xí)í = |0). 

The 2-qubit master equation now takes the following form: 

p = -i[H,p]+ J2 l/^ ) P^ ) ^p)+r{\+) l (+\t^p-p) 

t=l,2 

(2) 

For r = 0, the steady state of this master equation is diag- 
onal in the product basis |si)|s2), with si, S2 £ {0, 1}, and 

ai^\sj) = (— l) 3 '|sj) (z-basis), and there is no entanglement. 
The dephasing channel, described by the 7-part, eventually 
destroys all entanglement that may initially exist or may be 
created for a short period of time (small compared to 7" 1 ). 

For r — > 00 we can neglect the Hamiltonian and decoher- 
ence parts, and the density matrix will be quasi-permanently 
projected into the product of two |+)-states, and again there 
will be no entanglement between the qubits. The ques- 
tion is this: Between these two extremes — a pure prod- 
uct state on the one hand and a classically correlated state 
on the other hand — is there a combination of parameters 
g (Hamiltonian), 7 (decoherence), and r (reset) such that the 
steady state will be entangled? Now we show that this is in- 
deed the case. 

To solve the master equation (0, we expand 
the density operator in the z-basis: p(t) = 
J2 1 s ' 1 , s ' 2 , SuS2 =o C ^,s' 2 is 1 , S2 (t)\s' 1 ,s' 2 )(s 1 ,s 2 \. In this ba- 
sis, the master equation takes the form of the following 
16 coupled, linear, differential equations for the coeffi- 
cients: C s{ ^. SltS2 = {-ig[(-l) s '^ - (-1)^] + 
7 /2[(-l)«i+*i + (-1)4+^ _ 2] _ 2,-K .... + 

r /2{Co,^;o,s2 -f Ci )S ^ ; i )S2 + C s ' it o- iSl ,o + Cs'^ijsi.i}' Since 
the solutions are lengfhy we discuss here only the structure. 
The differential equations for the diagonal elements form 
a closed set. For any r > their value is 1/4 in steady 
state. For r = the steady state is not unique since the 
dephasing channel leaves the diagonal elements untouched, 
but the entanglement between the two qubits is zero for 
any initial vàlues. In the following, we will always assume 
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r > 0. The off-diagonal coefficients that are not on the 
anti-diagonal are coupled among themselves and are driven 
by the diagonal elements. Their final vàlues for t — > oo are 

n — C< — n* — n* — r(-ig+r+7 /2) 

^0001 — L·ooiO — L-Olll ~ °1011 — 4(2 9 2 + (r+ 7 /2)(r+ 7 ))- 

Finally, the anti-diagonal coefficients are not coupled among 
themselves and are driven by the other off-diagonal ele- 
ments. In steady state, the anti-diagonal is real and contains 

four times the following entry: ,, , w/iitV 2 !,,, , » . 

fe J 4( r + 7 )(2 ff ^ + (r+7/2)(r+7)J 

All other matrix elements are fixed by the Hermiticity 
of the density matrix. All coefficients approach steady 
state exponentially fast with their own characteristic 
exponents determined by the spectrum of the total Liou- 
ville super-operator C defined by p — Cp with vàlues: 
{0, -r, -2r, -2(r + 7 ),-(3r/2 + 7 ± 2i^Jg 2 - r 2 /16)} 
and multiplicities {1, 2, 1, 4, 4 + 4} respectively. 

As a measure for the entanglement between the two qubits 
we use the negativity @|, which is defined with respect to a 
bipartition A-À of the set of qubits as Ma = {\\p Ta I |i - l)/2, 
where Ta means the partial transpose with respect to A. For 
two qubits, the negativity can assume vàlues between zero 
(separable state) and 1/2 (maximally entangled state). 

In steady state, we compute from the above expressions for 
the density matrix the following analytic expression for the 
negativity in terms of the parameters g (Hamiltonian interac- 
tion), 7 (strength of the dephasing channel), and r (reset rate): 



7 (r + 7/2) 2 + g 2 (r + 7) - rg(r + 7) 

m X1 ' 2(r + 7)[2 5 2 + (r + 7/2)(r + 7)] * 

(3) 

Equation (|3} contains the full information about the entangle- 
ment properties of the two qubits and is one of the main re- 
sults of this paper. Note that M = M(g, r) depends only on 
two parameters, g = .9/7 and f = r/7. 

In Fig.^we see a contour plot of the negativity function. 
The key feature is the color-coded region in the r-<?-plane 
with steady state entanglement, where a darker color indicates 
higher entanglement. Three lines are marked in the entangled 
region: i) The upper white line is the maximum in g-direction 
(at constant r), ii) the lower white line is the maximum in r- 
direction (at constant g), and iii) the middle line is the straight 
g = r/(l + V3)- This middle line is approached by the upper 
and lower one asymptotically for large g,r. The global max- 
imum of the negativity is on this middle line at infinity with 
a value of approximately 0.092, about 20% of the maximal 
value. The darkest, most entangled area in our plot has neg- 
ativity 0.068. The entangled region is bounded by the grey 
line given by one of the roots of the non-trivial part of equa- 
tion (J3}- Outside of this region, the state is separable (white 
area). The entangled region approaches asymptotically the 
straights .9 = 7 and g = r plotted black in Fig.0 The asymp- 
totic line g = 7 is independent of r and indicates the border 
between the weak coupling and the strong coupling regime. 
For weak coupling, g < 7, decoherence/noise will always 
triumph over the Hamiltonian part that tries to create entan- 
glement. That is, as a necessary condition, we need to be in 
the strong coupling limit to observe entanglement. The inlet 
in Fig.^shows a cut at g = 2.57 through the color-plot. Most 



notable is the existence of a threshold value for r/7 above 
which entanglement is present in the steady state. 




FIG. 1: Separable states (white area) and entangled states (colored 
area) in the r-g-plane, where r is the rate of our reset process, and g is 
the coupling strength in the Ising Hamiltonian, and we use 7 -1 = 1 
as unit timescale. The color encodes the amount of entanglement 
measured by the negativity: the darker the area, the more entangle- 
ment is present. For a discussion of the other lines, please see the 
text. The inlet shows a cut for constant g — 2.57. 



We also simulated the system on a computer, where we use 
a microscopic decoherence model rather than a master equa- 
tion description. Decoherence is thereby induced by interac- 
tions of system partides with an environment, modelled as 
a background spin gas 01, and we consider interactions of 
Ising-type, which leads to dephasing noise. Since the spin gas 
exhibits memory, the decoherence process is non-Markovian 
and possibly non-local. Also in this case, the results qualita- 
tively agree with the analytic solution of our toy model. We 
remark that the system itself may also be described as a spin 
gas, even with fluctuating particle number. 

All these results support the following conclusion: Entan- 
glement can persist even in open noisy quantum mechanical 
systems that are far from thermodynamic equilibrium. It is 
true that even for r = 0, i.e. without reset mechanism, one 
can obtain entangled steady states, e.g. for a decay chan- 
nel (s = 1,-6/2 = C = 7) and an interaction Hamilto- 
nian H = gax &x ■ However, this phenomena is restricted 
to specific Hamiltonians and noise channels, where a steady 
state of the noise channel with sufficiently low temperature is 
a necessary condition for steady state entanglement. In partic- 
ular, entanglement vanishes for high temperatures of the bath, 
s — > 1/2. In contrast, we show in the next section that steady 
state entanglement appears generically in systems which fea- 
ture some kind of (appropriate) reset mechanism. 

Beyond toy model: many-particle systems. To demon- 
strate that we have identified a genèric feature, we will gen- 
eralize the system in various directions. Although the details 
will change, the qualitative behavior of the two-qubit model 
does not depend on the particular choice of the interaction 
Hamiltonian or the decoherence model. Fig. a) shows e.g. 
steady state entanglement for a XYZ-Hamiltonian as function 
of reset rate r, and decoherence described by a noise opera- 
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tor £ no i S e with different parameters. The qualitative behavior 
is similar to the inlet of Fig. [2 and we observe steady state 
entanglement even for infinite temperature of the bath. 

N (o) N (6) 



0.08 




r/7 r/7 

FIG. 2: (a) Negativity for 2-qubit system with XY Z-interaction and 
magnètic field, H = g^^a^ + OSa^a^ + aí 1} aí 2) + 
0.5(os + <J^)), as a function of the reset rate r at g = IO7. 
The noise is described by C n0 ise of equation Q with C = B/2 
and 7 = C/10. The upper curve corresponds to zero temperature 
(s = 0), the lower one to infinite temperature (s = 1/2) of the 
bath. Curves for any finite temperature lie in between. (b) Average 
negativity of (i) 5-qubits that have all pairwise interacted with each 
other (dashed), (ii) the reduced density matrix for any two qubits in 
a 5-qubit setting (dashed-dotted), and (iii) Poissonian mixture of re- 
duced density matrices corresponding to fluctuating particle number, 
see text, (sòlid) as a function of the reset rate r for Ising interaction 
and dephasing, with g — 57. 

Second, the idealized reset mechanism we consider can be 
replaced by a more realistic imperfect reset mechanism. In 
this case, fresh partides come in mixed states with sufficiently 
low entropy rather than in puré states (with entropy 0). Still, 
the steady state turns out to be entangled. 

Third, we have also considered systems consisting of more 
partides N, In this case we are interested in (i) multipar- 
ticle entanglement, and (ii) entanglement of reduced states 
of two qubits. In case of (i), we use the average negativity 
Sí, which is the negativity averaged over all possible biparti- 
tions of the system 0] to measure multiparticle entanglement, 
Non-zero average negativity ensures the presence of some 
form of entanglement in the system. To compute this quantity 
is a formidable task. The system size, and hence the number 
of differential equations and relevant bipartitions scales expo- 
nentially. To simplify computation, we consider a symmetric 
situation, where all qubits interact pairwise via Ising interac- 
tions, and compute Aí of up to seven qubits. Regarding (ii), 
we consider pairwise entanglement between two partides in 
7V-qubit systems, where the state of the two qubits is given by 
the reduced density matrix that is obtained by tracing out the 
remaining partides. Although typically such reduced-state 
entanglement is rather unstable in multiparticle entangled sys- 
tems (since entanglement with other system partides acts as 
additional noise source), we find even in this case a parameter 
regime with entanglement in steady state. However the re- 
gion in r — g diagram where one finds entanglement becomes 
smaller for increasing number of partides N. Fig.|2jb) shows 



the plot of the average negativity and negativity of reduced 
states for systems of N = 5 qubits. 

Finally, we have considered systems with a fluctuating 
number of partides N. To be precise, we consider sys- 
tems of N qubits that interact pairwise in a symmetric way, 
and TV fluctuates between 2 < N < 6 following a Pois- 
sonian distribution in this interval. The state of two qubits, 

Pab = J2n PnPab-* i s gi yen by a mixture of reduced density 
operators p^2 (obtained from the steady state of the ./V-qubit 
system) with mixing probabilities given by the probability 
that the system size is N. We find that entanglement between 
pairs of qubits remains finite (see Fig.|2jb)). 

Generally, for any interacting system subjected to some 
kind of noise or decoherence (not necessarily local or Marko- 
vian), one can expect that an appropriate local reset mech- 
anism leads to entangled steady states, as long as the sys- 
tem is capable of creating entanglement from (slightly mixed) 
input states at some point in time. For reset rates that are 
slow enough that entanglement may be built up in the system, 
but fast enough that decoherence (or other noise processes) 
will not destroy all the entanglement, one typically expects 
steady state entanglement on an appropriate coarse grained 
time scale. 

Note that various physical processes are conceivable for the 
reset mechanisms, including active processes such a measure- 
ment or "optical pumping" of randomly selected partides in 
a spin gas, as well as simply replacing a qubit by a fresh 
one (e.g. spin dependent tunnelling of electrons in charge 
controlled quantum dots). A process of the latter type may 
also be conceivable in more natural scenarios such as bio- 
molecular interactions in a celi. These systems are far away 
from the thermodynamical equilibrium, and there the fluctua- 
tion of partides may, at the same time, act as reset mechanism. 

Summary: We have shown that entanglement can be 
present in open noisy quantum systems far from thermody- 
namic equilibrium. For a two-qubit toy model, we could an- 
alytically solve the master equation consisting of a Hamilto- 
nian part, a noise channel, and a proposed reset mechanism, 
which is responsible for the non-vanishing entanglement in 
the steady state. We were able to give a closed expression 
for the entanglement as a function of the parameters of the 
master equation. We extended the analysis to similar models 
with other interaction Hamiltonians and decoherence models 
(including non-Markovian decoherence), systems consisting 
of more partides, and also systems with fluctuating particle 
number. In all cases we found that steady state entanglement 
can prevail. This demonstrates the possibility of stable steady 
state entanglement in natural systems consisting of many par- 
tides which are far away from thermodynamical equilibrium. 
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